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PROGRAM  DESCRIPTION 


INTRODUCTION 


The  Solid  Propulsion  Optimization  Code  (SPOC)  performs  detailed 
preliminary  designs  of  a  large  variety  of  solid  propellant  rocket  motors. 
Dimensions  of  the  propellant  grain,  nozzle,  and  pressure  vessel  are 
adjusted  by  the  code,  along  with  propellant  formulation  and  burn  rate,  to 
produce  a  motor  design  that  satisfies  performance  requirements  and 
design  constraints,  and  that  has  been  optimized  with  respect  to  a  user- 
selected  parameter. 

This  volume  of  the  User's  Manual  -  Volume  III  (Program  Description)  - 
contains  the  subroutine  descriptions  and  flow  charts,  and  cross-indices  of 
common  statements,  subroutines  and  call  statements.  Volume  I  (Technical 
Description)  gives  the  basis  for  the  code  computations  analytical  develop¬ 
ment,  logic  flow  charts  used  in  verification  checks  and  error  messages. 
Volume  II  (User's  Guide)  contains  the  input  and  output  dictionaries  and  their 
accompanying  illustrations,  along  with  other  input  instructions  needed  to 
execute  the  code. 
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SUBROUTINE  DESCRIPTIONS 


ABMODL 

ACCEL 

ACOUST 

AF2SUB 

ARC 

AREAS 

ASESUB 

AS2SUB 

AS3SUB 

ATAESB 

ATMOS 

AWESUB 


Calculates  combustion  response  per  analytical  response 
model.  Used  by  subroutine  RSPNSE. 

Calculates  missile  horizontal  and  vertical  acceleration.  A 
point-mass  missile  is  assumed  to  be  flying  a  ballistic  tra¬ 
jectory  in  the  altitude -range  plane. 

Calculates  velocity  of  acoustic  volume  current  at  end  of 
motor.  Used  by  subroutine  MODCLC. 

Solves  for  port  perimeter  and  area  of  Type  2  configurations 
for  each  plane.  This  is  a  2-dimensional  solution  and  is 
called  by  LP2SUB.  See  Note  1. 

Surface  area  computational  subroutine  used  with  forward 
propellant  segment  of  conocyl  grain. 

Surface  area  computational  subroutine  used  with  forward 
propellant  segment  of  conocyl  grain.  Calls  LOOK1  and 
LOOK2,  which  are  separate  entries  in  subroutine  ELIPS. 

Calculates  surface  areas  and  volumes  of  individual  segments 
and  end  sections  by  calling  RASUB,  XRSUB,  LRSUB  and 
HESUB.  These  break  the  configuration  down  into  equivalent 
dimensions  for  a  simplified  configuration. 

First  level  geometry  for  end  closures  of  type  2  configura¬ 
tions.  This  is  a  3-dimensional  calculation.  Being  primarily 
an  executive  routine,  it  calls  ASESUB,  AWESUB  and  RCSUB. 
See  Note  1. 

First  level  geometry  for  end  closures  of  type  3  configura¬ 
tions.  This  is  a  3- dimensional  calculation.  Being  primarily 
an  executive  routine,  it  calls  ASESUB,  AWESUB,  HESUB, 
and  RCSUB.  See  Note  1. 

Calculates  nozzle  throat  area  (AT)  and  exit  area  (AE  at 
any  time  during  ballistic  simulation. 

Calculates  atmospheric  properties  as  a  function  of  altitude 
for  (1)  1959  ARDC  STD  DAY;  (2)  MIL-STD-210A  TROPICAL 
DAY;  (3)  MIL- STD- 2 10 A  POLAR  DAY;  (4)  MIL-STD-210A 
HOT  DAY;  (5)  MIL-STD-210A  COLD  DAY. 

Calculates  areas  and  volumes  for  the  web  segment  of  a 
type  2  configuration.  A  3- dimensional  end  routine. 

See  Note  1. 
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Subroutines  in  SPOC  (Continued) 


BLDATA 

CASEAN 

CFVLSB 

CHEKIN 

CLOS 

COMP 

CPHS 

CROSS1 

CROSS2 

D  PRINT 

ECSUB 

ELIPS 

EQLBRM 

EXIT 


Biock  data  for  thermochemistry:  Atomic  Symbols,  Atomic 
Weights,  Valences. 

Calculates  case  wall  thickness,  forward  and  aft  closure 
thicknesses  and  weight  summary  of  case  structural  parts. 

Calls  WACL1,  WFCL1,  and  WF2A3,  which  are  separate 
entries  in  subroutine  CLOS. 

Calculates  Cpv^  (thrust  coefficient  for  vacuum  conditions 
corrected  for  divergence  losses).  Xn  is  brought  into 
■his  subroutine. 

An  input  checking  and  diagnosis  routine.  Prints  messages 
that  describes  the  problem  the  user  has  Bet  up. 

Case  forward  and  aft  closure  computational  routine. 

Calculates  geometric  details  for  three  forward  closure  types 
and  two  aft  closure  types,  plus  associated  volumes  and  weights. 

Executive  (Main  Logic  Control)  routine.  Calls  all  other 
analysis  modules. 

Calculation  of  thermodynamic  properties  of  individual 
species.  Part  of  thermochemistry  module. 

Sets  up  for  subroutine  ACOUST.  Used  by  subroutine  MODCLC. 

Detects  sign  change  (zero)  of  acoustic  volume  current 
velocity  to  within  convergence  tolerances.  Used  by  subroutine 
MODCLC. 

D-array  summary  print  routine.  The  D-array  is  the  list  of 
variables  upon  which  the  optimizer  is  allowed  to  operate. 

Calculates  and  tests  constants  that  deal  with  the  end  closure; 
primarily  tests. 

Grain  regression  computational  routine  for  conocyl  grain 
with  elliptical  forward  closure. 

Calculates  equilibrium  composition  and  properties  of  propel¬ 
lant  combustion  species.  Part  of  thermochemistry  module. 

Universal  exit  to  give  print-out  of  current  values  of  PATSi.1 
adjusted  parameters  when  there  is  an  abnormal  termination 
of  code  execution. 
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Subroutines  in  SPOC  (Continued) 


E488M2 

FIGS 

FIG6 

FIG7 

FLT 

FPACC 

GAUSS 

G2F2FG 

HESUB 

HITEMP 

IB  SUB 

IMPEFF 

INPUT 

INTGE 

ITERP1 


Executive  subroutine  for  stability  analyses.  Calls  all  other 
subroutines.  Calculates  and  writes  stability  penalty. 

Table  look-up  for  photoelastic  stress  parameter  HSTAR  ver¬ 
sus  valley  included  angle  to  UBe  with  star  grain  configuration. 

Table  look-up  for  photoelastic  stress  parameter  HSTAR  ver¬ 
sus  valley  width  to  use  with  star  grain  configuration. 

Table  look-up  for  photoelastic  stress  parameter  HSTAR  ver¬ 
sus  valley  included  angle  to  use  with  wagon-wheel  grain  con¬ 
figuration. 

Calculates  ideal  drag -free  velocity  at  end  of  motor  operation 
and  maximum  axial  acceleration  during  motor  operation. 

Part  of  missile  flight  dynamics  logic.  A  diagnostic  routine 
used  in  the  computation  of  missile  acceleration  along  the 
flight  path. 

Solution  of  N  simultaneous  linear  equations.  Used  in  ther¬ 
mochemistry  module. 

Calculates  functions  used  in  calculation  of  stability  integrals. 
Used  by  subroutine  INTGL. 

Calculates  the  distance  from  the  last  cylindrical  portion 
plane  to  the  end  of  the  end  section  for  given  distances  from 
motor  centerline. 

Computes  motor  operating  limit  parameters  with  high  tem¬ 
perature  grain  and  compares  them  to  requirements. 

Internal  Ballistics  Subroutine  -  Calculates  all  gas  flow  pro¬ 
perties  between  any  two  adjacent  planes.  This  is  the  heart 
of  subroutine  SEC3SB. 

Computes  rocket  motor  impulse  efficiency  using  SPP 
empiricisms  (Ref.  AFRPL-TR-75-36). 

Reads  namelist  STABIN.  Sorts  incoming  data  set  from 
common/TIMDAT/.  Calculates  pressure  and  velocity 
exponents.  Prints  input  data. 

Calculates  stability  integrals.  Used  by  subroutine  E488M2. 

Nth  order  interpolation  of  a  table  having  a  single  independent 
variable. 
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Subroutines  in  SPOC  (Continued) 


LIQUID 

LOTEMP 

LP1SUB 

LP2SUB 

LP3SUB 

LP4SUB 

LRSUB 

MACHNO 

MAIN  CO 

MATRIX 

MODCLC 

MTRCST 

NEWING 

NEWRAP 


Part  of  thermochemistry  input  diagnostics.  Checks  numeri¬ 
cal  compatibility  of  amounts  of  liquid  ingredients  inpui  and 
total  solids  level. 

Computes  motor  operating  limit  parameters  with  low  tem¬ 
perature  grain  and  compares  them  to  requirements. 

Calculates  the  port  perimeter  and  area  for  type  1  configura¬ 
tions  at  each  plane.  See  Note  1. 

Calculates  the  port  perimeter  and  area  for  type  2  configura¬ 
tions  at  each  plane.  See  Note  1. 

Calculates  the  port  perimeter  and  area  for  type  3  configura¬ 
tions  at  each  plane.  See  Notte  1. 

Calculates  the  port  perimeter  and  area  for  type  4  configura¬ 
tions  at  each  plane.  See  Note  1. 

Calculates  the  section  of  the  perimeter  between  radii  in 
the  end  section. 

Solution  for  mach  number  from  area  ratio  and  specific  heat 
ratio.  Ideal,  one-dimensional,  single  phase  gas  dynamics 
are  assumed. 

Executive  {main  control)  routine  for  thermochemistry  compu¬ 
tation  once  the  incoming  ingredient  weight  fraction  set  has 
been  checked  for  numerical  validity  and  for  compatibility 
with  user-input  limits. 

Matric  solution  used  in  thermochemistry  module. 

Detects  acoustic  mode  frequency  by  convergence  on  zeros 
of  acoustic  volume  current  velocity.  Used  by  E488M2. 

Calculates  development,  qualification,  PFRT  and  production 
costs  for  steel  case  motors  using  Tri-Services  cost  model. 

Manipulative  routine  for  modification  of  thermochemistry 
logic  by  adding  new  ingredients 

Newton- Raphson  iterative  solution  of  an  implicit  function. 
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Subroutines  in  SPOC  (Continued) 


NORMAL 

Adjusts  incoming  set  of  propellant  ingredient  weight 
fractions  so  that  total  equals  1.  0.  Adjusts  only  those 
ingredients  that  user  has  specified  for  adjustment  during 
optimization  process. 

NOZINP 

Nozzle  configuration  data  reading,  check,  and  print. 
Determines  geometric  validity  of  incoming  data  set. 

NOZL 

Calculates  nozzle  geometry  from  incoming  data  set. 

Performs  thermal  and  structural  analyses.  Calculates 
weight. 

NOZZSB 

Lists  nozzle  parameters  used  by  ballistic  simulation  module 
and  checks  them.  It  is  called  by  subroutine  SEC1SB. 

NREQ1 

Iterative  solution  for  a  dimension  (ALPHMX)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ2 

Iterative  solution  for  a  dimension  (AFSTAR)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ3 

Iterative  solution  for  a  dimension  (AFSTAR)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ4 

Iterative  solution  for  a  dimension  (BETAX)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ5 

Iterative  solution  for  a  dimension  (DELTA)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ6 

Iterative  solution  for  a  dimension  (GAMMA)  of  Type  1  (Star) 
grain  required  by  ballistic  simulation  module. 

NREQ7 

Iterative  solution  for  a  dimension  (GAMMA)  of  Type  2  (Wagon 
Wheel)  grain  required  by  ballistic  simulation  module. 

NREQ8 

Iterative  solution  for  a  dimension  (GAMMA)  of  Type  2  (Wagon 
Wheel)  grain  required  by  ballistic  simulation  module. 

NREQ9 

Iterative  solution  for  a  dimension  (GAMMA)  of  Type  2  (Wagon 
Wheel)  grain  required  by  ballistic  simulation  module. 
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Subroutines  in  SPOC  (Continued) 


ONETMP 

OUT1 

PATSH 

PCGSUB 

PC1SUB 

PC2SUB 

PC3SUB 

PC4SUB 

PEPCSB 

PROPST 

PTLSTB 

RASUB 


Calculates  operating  limits  and  associated  penalties  when 
ballistic  simulation  is  performed  at  only  one  grain  temper¬ 
ature. 

Converts  part  of  thermochemical  calculation  results  to 
proper  units  and  rocket  definitions.  Called  by  RKTOUT. 

Pattern  search  optimization  routine. 

Is  a  general  solution  of  plane  constants  for  either  side  of 
a  Type  2  configuration.  It  does  the  calculations  for  PC2SUB. 

See  Note  1. 

Prints  input  table  of  "Thickness,  Perimeter  and  Port  Area"; 
checks  geometry  values  and  completes  geometry  print-out; 
calls  LP1SUB;  and  prints  "Initial  Perimeter,  Initial  Port 
Area  and  Inert  Sliver"  values  for  Type  1  configurations. 

All  of  this  is  seen  in  geometry  part  of  the  print-out.  See  Note  1. 

Prints  out  the  geometry  section  of  the  input  data  for  Type 

2  configurations.  This  routine  is  primarily  executive.  It 
calls  LP2SUB  and  PCGSUB.  See  Note  1. 

Prints  out  the  geometry  section  of  the  input  data  for  Type 

3  configurations.  See  Note  1. 

Prints  out  geometry  input  data;  calls  LP4SUB  and  prints 
out  "Initial  Perimter,  Initial  Port  Area  and  Inert  Sliver" 
data  for  Type  4  configurations.  See  Note  1. 

Calculates  ratio  of  pressure  at  the  nozzle  exit  to  pressure 
in  the  chamber.  This  is  an  iterative  solution. 

Performs  grain  structural  analysis. 

Calculates  stability  margin  (linear  alpha)  and  damping 
coefficient.  Prints  stability  summary  and  details.  Used 
by  subroutine  E488M2. 

Calculates  the  radius  to  any  point  on  the  burning  surface 
perimetdr  in  the  end  sections. 
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Subroutines  in  SPOC  (continued) 


RATES B 

RCSUB 

REACT 

RKTOUT 

ROCKET 

RSPNSE 

RTCKSB 

SAVE 

SBFHAT 

SBQBAR 

SEARCH 

SEC1SB 


Calculates  burn  rate  as  a  function  of  pressure  and  Mach 
number.  Called  mainly  by  subroutine  IBSUB. 

Calculates  radius  from  the  motor  centerline  to  the  closure 
in  the  end  sections. 

Locates  propellant  ingredients  in  internal  tables  in  thermo¬ 
chemistry  module  and  sets  proper  indices  for  locating 
elements  and  coefficients. 

Converts  part  of  thermo  chemical  calculation  results  to  proper 
units  and  rocket  definitions.  Calls  OUT3,  a  second  entry  to 
subroutine  OUT1. 

Converts  part  of  thermochemical  calculation  results  to  proper 
units  and  rocket  definitions.  Calls  NEWOF,  a  second  entry 
to  subroutine  SAVE. 

Calculates  combustion  response  per  selected  response  option 
(four  different  models)  for  each  section  of  the  motor.  Used 
by  subroutine  E488M2. 


Checks  for  validity  of  the  burn  rate  inputs. 

Saves  thermochemical  data  for  subsequent  passes  with  same 
ingredients  to  speed  solution. 

Calculates  the  acoustic  pressure  at  a  specific  location  in 
the  motor.  Used  by  subroutine  STBINT. 

Calculates  the  mean  flow  volume  current  at  a  specific 
location  in  the  motor.  Used  by  subroutine  STBINT. 

Looks  for  thermodynamic  coefficients  of  combustion 
species  using  indices  created  by  subroutine  REACT. 

Determines  the  number  of  input  planes.  Interpolates 
between  plane  inputs.  Computes  plane  constants  by  calling 
PC(1,2,3  or  4)  SUB  depending  on  the  type  geometry.  Com¬ 
putes  head  end  and  nozzle  end  constants  by  calling  ECSUB. 
Computes  initial  volumes.  Checks  scale  factor  inputs  and 
sets  SF  =  1,  0  if  scale  factors  are  not  used.  Checks  nozzle 
diameter  inputs  by  calling  NOZZSB.  Checks  time  incre¬ 
ment  inputs.  Checks  burning  rate  inputs  by  calling  RTCKSB. 
Prints  all  the  values  it  calculated  or  checked. 
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Subroutines  in  SPOC  (Continued) 


SEC2SB 


SEC3SB 


SETUP1 


SETUP2 


SETUP3 


SETUP4 


SETUP5 


STBDAT 


Computes  end  section  surface  areas  and  volumes  by  calling 
LP(2  or  3 )SUB  and  AS(2  or  3)SUB  depending  on  the  geometry. 
Computes  perimeters  and  port  areas  by  calling  LP(1,2,3  or 
4)SUB  depending  on  the  geometry.  Calculates  fuel  areas, 
incremental  surface  areas,  and  propellant  volumes.  Com¬ 
putes  interpolations  on  burn  rate  scale  factors  if  they  are 
used.  Computes  nozzle  throat  and  exit  areas  by  calling 
ATAESB.  Computes  ballistic  performance  by  calling 
SEC3SB.  Updates  time  and  web  thickness  for  next  incre¬ 
ment.  Goes  back  to  the  first  of  the  subroutine  for  the 
next-time  calculations. 

Calculates  ballistic  constants.  Computes  ballistics  for 
the  head-end  and  nozzle-end  segments  by  calling  RATESB 
and  IBSUB  in  sets.  Cojnputes  nozzle  section  ballistics 
as  in  the  head  end.  Checks  mass  generated  with  mass 
discharged  for  equilibrium.  Adjusts  pressure  for  new 
run  if  equilibrium  is  not  reached.  Computes  thrust  and 
integrals.  Resets  initial  conditions  for  the  next  time 
increment.  Statement  50  is  the  beginning  of  the  main 
ballistic  loop. 

Input  validity  tests  and  geometry  setup  for  communication 
with  ballistic  simulation  module  for  grain  Type  1  (Star). 

See  Note  2. 

Input  validity  tests  and  geometry  setup  for  communication 
with  ballistic  simulation  module  for  grain  Type  2  (Wagon 
Wheel).  See  Note  2. 

Input  validity  tests  and  geometry  setup  for  communication 
with  ballistic  simulation  module  for  grain  Type  3  (Finocyl). 
See  Note  2. 

Input  validity  tests  and  geometry  setup  for  communication 
with  ballistic  simulation  module  for  grain  Type  4  (Conocyl). 
See  Note  2. 

Input  validity  tests  and  geometry  setup  for  communication 
with  ballistic  simulation  module  for  grain  Type  5  (CP). 

See  Note  2. 

Creates  matrix  of  grain  dimensions  and  internal  gas  dynamic 
conditions  for  use  by  combustion  stability  module.  Called  by 
SEC3SB.  Calls  STBCLO  and  STBNOZ  (which  are  separate 
entries  into  CLOS  and  NOZL,  respectively). 
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Subroutines  in  SPOC  (Contined) 


STBINT 

TCHEM 

TRAJ 

TRAJIN 

UBCALC 

USERCS 


Calculates  the  stability  integrals  for  each  section  of  the 
motor.  Used  by  subroutine  E488M2. 

Executive  subroutine  for  propellant  thermochemical 
analysis  module.  Called  by  subroutine  COMP. 

Trajectory  flight  dynamics  computational  routine.  Point 
mass  missile  flying  two  dimensional  ballistic  trajectory  in 
altitude  range  plane.  Fourth  order  Runge-Kutta  numerical 
integration. 

Trajectory  data  read,  decision  logic,  and  print  routine. 

Calculates  constants  describing  geometry  (perimeter  and 
port  area)  and  mean  flow  volume  current  for  each  section 
of  the  motor.  Used  by  subroutine  E488M2. 

User-supplied  cost  model. 


USEREF  User- supplied  impulse  efficiency  model. 


USERRB  User- supplied  propellant  burn  rate  model. 

USERRH  User- supplied  propellant  rheological  properties  model. 

USERSE  User- supplied  propellant  strain  endurance  model. 


VOLUME  Component  volume  computational  routine  for  nozzle 

analysis  subroutine  NOZL. 

XRSUB  This  routine  calculates  the  distance  (XR)  from  the  boundary 

of  a  symmetrical  segment  to  the  perimeter.  This  is  used 
by  the  end  section  subroutines. 


NOTES 

1.  Grain  type  numbers  in  noted  definitions  are  designations  internal  to  the 
ballistic  simulation  module  and  are  not  the  same  as  the  grain  types 
available  in  SPOC. 

2.  Calls  AFTCL1,  AFTCL2,  FCL2A3,  and  FWDCLl,  which  are  separate 
entries  in  subroutine  CLOS. 
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SUBROUTINE  CROSS-INDEX 


Subroutine 

Called 

by  These 

Subroutines 

ABMODL 

RSPNSE 

ACCEL 

TRAJ 

ACOUST 

MODCLC 

AFTCL1 

E 

CLOS 

SETUPl 

SETUP2 

SETUP3 

SETUP4 

AFTCL2 

E 

CLOS 

SETUPl 

SETUP2 

SETUP 3 

SETUP4 

AF2SUB 

LP2SUB 

ARC 

AREAS 

AREAS 

SETUP4 

ASESUB 

AS2SUB 

AS3SUB 

AS 2 SUB 

SEC1SB 

SEC2SB 

AS 3 SUB 

SEC1SB 

SEC 2 SB 

ATAESB 

SEC2SB 

ATMOS 

ACCEL 

TRAJ IN 

AWE SUB 

AS 2 SUB 

AS 3 SUB 

BLDATA 

CASEAN 

COMP 

CFVLSB 

SEC3SB 

CHECKIN 

MAIN 

CLOS 

SETUPl 

SETUP 2 

SETUP 3 

SETUP4 

SETUPS 

COMP 

MAIN 

CPHS 

EQLBRM 

CR0SS1 

MODCLC 

CROSS2 

MODCLC 

DPRINT 

COMP 

ECSUB 

SEC1SB 

EL  I  PS 

AREAS 

EQLBRM 

ROCKET 

EXIT 

AWESUB 

CNTRL 

COMP 

ELIPS 

E488M2 

INPUT 

MODCLC 

NOZINP 

SEC3SB 

SETUPl 

SETUP 3 

SETUP4 

SETUP 5 

TCHEM 

TRAJIN 

E488M2 

COMP 

FCL2A3 

E 

CLOS 

SETUPl 

SETUP2 

SETUP 3 

SETUP5 

FIG5 

PROPST 

FIG6 

PROPST 

FIG7 

PROPST 

FLT 

COMP 

FPACC 

TRAJ 

FWDCLl 

E 

CLOS 

SETUPl 

SETUP2 

SETUP3 

SETUP 4 

GAUSS 

EQLBRM 

G2F2FG 

INTGL 

HESUB 

ASESUB 

AS3SUB 

AWESUB 

HI TEMP 

COMP 

IBSUB 

SEC3SB 

IMPEFF 

COMP 

INPUT 

E488M2 

INTGL 

E488M2 

* 

ITERP1 

ACCEL 

EITEMP 

IMPEFF 

ITERP2 

LOTEMP 

RSPNSE 

LIQUID 

TCHEM 

LOOK  2 

E 

AREAS 

EL  I  PS 

LOOK  3 

fi* 

AREAS 

ELIPS 
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SETUPS 

SETUPS 


FPACC 

SETUP2 


SETUP5 


ONETEMP 


Subroutines  Cross  Index  (Continued) 
Subroutine  Called  by  These  Subroutines 


LOTEMP 

COMP 

LP1SUB 

PC1SUB  SEC2SB 

LP2SUB 

PC2SUB  SEC1SB 

LP3SUB 

PC3SB  SEC1SB 

LP4SUB 

PC4SB  SEC 2 SB 

LRSUB 

ASESUB 

MACHNO 

NOZL 

MAIN 

MAINCO 

TCHEM 

MATRIX 

EQLBRM 

MODCLC 

E488M2 

MTRCST 

COMP 

NEWING 

TCHEM 

NEWOF  E 

ROCKET  SAVE 

NEWRAP 

SETUP1 

NORMAL 

TCHEM 

NOZINP 

COMP 

NOZL 

COMP 

NOZZSB 

SEC1SB 

NREQl 

SETUP 1 

NREQ2 

SETUP 1 

NREQ3 

SETUP 1 

NREQ4 

SETUP 1 

NREQ5 

SETUP 1 

NREQ6 

SETUPl 

ONETEMP 

COMP 

OUT1 

RKTOUT 

OUT  3  E 

OUT1  RKTOUT 

PATSH 

MAIN 

PCGSUB 

PC 2 SUB 

PC1SUB 

SEC1SB 

PC2SUB 

SEC1SB 

PC3SUB 

SEC1SB 

PC4SUB 

SEC1SB 

PEPCSB 

SEC 3 SB 

PROPST 

COMP 

PTLSTB 

E488M2 

RASUB 

ASESUB 

RATESB 

HI TEMP  IB SUB 

RATESB 

INPUT 

RCSUB 

AS 2 SUB  AS3SUB 

REACT 

MAINCO 

RKTOUT 

ROCKET 

ROCKET 

MAINCO 

ROCKETl  E 

ROCKET 

RSPNSE 

E488M2 

RTCKSB 

SEC1SB 

SAVE 

ROCKET 

SBPHAT 

STBINT 

SBQBAR 

STBINT 

SEARCH 

MAINCO 

SEC 1 SB 

COMP 

SEC2SB 

SEC2SB 


LOTEMP 


ONETMP 


SEC3SB 
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Subroutine  Cross-Index  (Continued) 


Subroutine  Called  by  These  Subroutines 


SEC 2 SB 

COMP 

SEC3SB 

SEC 2 SB 

SETUPl 

COMP 

SETUP2 

COMP 

SETUP3 

COMP 

SETUP4 

COMP 

SETUPS 

COMP 

STBCLO 

CLOS 

STBDAT 

SEC3SB 

STBINT 

E488M2 

STBNOZ 

NOZL 

STBS2 

STBDAT 

TCHEM 

COMP 

TRAJ 

COMP 

TRAJIN 

COMP 

UBCALC 

E48SM2 

USERCS 

COMP 

USEREF 

COMP 

USERRB 

HI TEMP 

USERRH 

TCHEM 

USERSE 

PROPST 

VOLUME 

NOZL 

WACLl 

E 

CASEAN 

WFCLl 

E 

CASEAN 

WF2A3 

E 

CASEAN 

XRSUB 

ASESUB 

STBDAT 

STBDAT 

TRAJIN 

IBSUB  LOTEMP 

CLOS 

CLOS 

CLOS 


ONETMP  SEC 3 SB 


NOTES 

(1)  Supplementary  notation  "E"  (e.g.  "AFTCL1  E")  indicates  an 
alternate  entry  into  the  subroutine  in  the  second  column. 
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Figure  1.  Subroutine  -  Common  Statement  Cross  Index  (Cont'd) 
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MAIN:  Reads  control  inputs;  initializes  some  parameters;  controls 

printout;  calls  search  routine 

COMP:  Executive  subroutine  passes  information  between  subroutines; 

calculates  some  penalties  and  overall  objective  function  (OBJ); 
provides  printout 

PATSH:  Adjusts  specified  parameters;  evaluates  changes  in  objective 

function  (OBJ) 


Figure  2.  Overall  Code  Organization 
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Figure  3.  Subroutine  COMP  Flow  Chart 
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Figure  4.  Subroutine  SETUPi  Flow  Chart 
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Figure  8.  Subroutine  E488M2  Flow  Chart 
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